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We propose a simple way to improve the laser gravitational-wave detectors sensitivity by means 
of reduction of the number of reflective coating layers of the core optics mirrors. This effects in the 
proportional decrease of the coating thermal noise, the most notorious among the interferometers 
technical noise sources. The price for this is the increased quantum noise, as well as high require- 
ments for the pump laser power and power at the beamsplitter. However, as far as these processes 
depend differently on the coating thickness, we demonstrate that a certain trade-off is possible, yield- 
ing a 20-30% gain (for diverse gravitational wave signal types and interferometer configurations), 
providing that feasible values of laser power and power on the beamsplitter are assumed. 



I. INTRODUCTION 

Among the multitude of noise sources limiting the sen- 
sitivity of contemporary laser gravitational- wave (GW) 
detectors (LIGO [1, 2], VIRGO [3, 4], GEO-600 [5, 6], and 
TAMA [7, 8]), the one, usually referred to as quantum 
noise and stemming from the quantum nature of light, 
stands apart from the rest of the noise sources, referred 
to as technical or classical ones, respectively. Quan- 
tum noise originates from quantum fluctuations of phase 
(shot noise) and intensity (radiation-pressure noise) of 
light circulating inside the interferometers, which obey 
Heisenberg's uncertainty relation [9] and therefore can 
not be reduced simultaneously. The latter group com- 
prises various fluctuations of thermal, seismic and sim- 
ilar origin that can be, in principle, diminished either 
by cooling, or using better materials, more sophisticated 
seismic isolation and so on. All the hitherto undertaken 
efforts towards the improvement of GW interferometers 
sensitivity went in two parallel but virtually independent 
streams. In many proposed methods of diverse techni- 
cal noise sources mitigation [10-18], the authors assumed 
quantum noise of the interferometer as independent of 
the technical noise budget and thus did not take it into 
consideration. Equally lukewarm were the researchers of 
the quantum noise, who proposed a plenty of sophisti- 
cated and witty ways for its reduction [19-33], towards 
the technical noise sources, tacitly implying them being 
independent on the quantum fluctuations of light inside 
the interferometer. 

However, the growing interest to the optimal config- 
urations of future GW detectors, inspired by the re- 
cent achievements in reduction of technical noise in first- 
generation detectors and a proximity of the start of con- 
struction of the second generation ones, has brought the 
problem of simultaneous treatment of the two groups of 
noise sources to the fore [29, 30]. Indeed, the sensitivity 
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of the current first generation detectors (the most sensi- 
tive of them, LIGO, has already finished its life cycle) is 
limited mostly by seismic noise at lower frequencies (be- 
low ~ 100 Hz) and by the quantum shot noise at higher 
frequencies. In the next-generation detectors, such as 
Advanced LIGO [34, 35], Advanced VIRGO [36], and 
LCGT [37], the technical noise will be reduced signifi- 
cantly by using much better seismic isolation and other 
technological advances. Along with this, the quantum 
shot noise will be suppressed by about one order of mag- 
nitude due to increased optical power and, very proba- 
bly, by the injection of the quantum squeezed light into 
the interferometer, as it has been proposed first in [9] 
(this technology has been successfully tested experimen- 
tally in GEO-600 [38]). Yet the second generation detec- 
tors sensitivity remain bound by the mix of the technical 
and quantum noise: shot noise at high frequencies, (ii) 
radiation-pressure noise at low frequencies, and (iii) coat- 
ings thermal noise in the best sensitivity medium frequen- 
cies band around 100 Hz (see Fig. 3 of the paper [35]). 

In this article, we make a further step forward and 
consider quantum noise and coating thermal noise in con- 
junction. The importance of thermal fluctuations in core 
optics dielectric coatings was realised by the community 
several years ago [39] . All the hitherto proposed ways to 
reduce these fluctuations, were it based on broadening 
of laser beams [10, 11, 14, 16], or on using better coat- 
ing materials [12], or on coating structure optimisation 
[13, 15, 17], followed implicitly a common rule that the 
end mirrors (ETMs) of the interferometer should be as 
reflective as possible thus requiring the number of coating 
layers to be pretty large (ca. 40). But since the power 
spectral density of these fluctuation rises linearly with 
the number of coating layers, the improvement provided 
by these methods is rather modest. 

Another apparent way to get rid of coating noise by 
getting rid of the (at least a part of) the coating itself 
was proposed even earlier. In [40], authors suggested to 
replace the end mirrors by coatingless corner reflectors, 
while in [41] the short anti-resonance-tuned Fabry-Perot 
cavities have to play the role of ETMs. Even more radical 
solution to use the pass-through Mach-Zehnder /Fabry- 
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Perot topology instead of the Michelson/Fabry-Perot one 
was suggested in [42]. However, the implementation of 
these methods in the near future is improbable, for corner 
reflectors were shown to have high optical losses [43] and 
the latter two solutions require too radical modiflcations 
of the GW detector optical setup. 

The requirement for the ETMs of the standard 
Michelson/Fabry-Perot scheme of GW interferometer 
(see Fig. 1) to have high reflectivity has rather strong 
logic behind. There are two obvious reasons in favour 
of it. First, in the power-recycled topology, the value of 
ETM power reflectivity i^ETM deflnes how much circu- 
lating power Ic can be built up in the arm cavities, for a 
given value of the input laser power /q: 

^ ~ 2(1 — R^tm) ^ ^ 

Simple estimate based on the values of /q and and Ic 
planned for the Advanced LIGO yields 1 — i^ETM ^ 10~^. 
Second, a non-ideal reflectivity of the end mirrors means 
an injection of additional optical vacuum fluctuations 
into the arm cavities, that is, the increase of the quantum 
noise. 

But if we ask ourselves a question whether this re- 
quirement always provides unconditionally optimal sen- 
sitivity for the GW interferometer, we claim that the 
answer will be 'no'. The circulating power issue can be 
solved either by using more powerful laser, as e.g. the 
one with /q = 500 W proposed for the third generation 
Einstein Telescope gravitation- wave detector [44], or by 
using squeezed vacuum injection as it might be done in 
Advanced LIGO [45]. 

Regarding the influence of the additional vacuum en- 
tering the arm cavities through the more transparent 
ETMs, it is evident, that in the scenario where the total 
noise budget, in the best sensitivity frequency band, is 
dominated by the mirrors coatings thermal noise, it is 
reasonable to reduce the number of coating layers of the 
end mirrors, increasing the quantum noise, but decreas- 
ing the coatings thermal noise. 

Therefore, the number of layers of the core optics coat- 
ings should be included into the set of optical parame- 
ters, such as arm cavities bandwidth, the signal recycling 
mirror transmittance and the the signal recycling cavity 
detuning, over which the minimisation of the sum noise 
of the interferometer is run. 

In this paper, we perform this kind of optimisation for 
the signal- and power-recycled Fabry-Perot-Michelson 
interferometer (see Fig. 1), assuming its main parame- 
ters close to the ones planned for the Advanced LIGO. 
In the next section, we describe the variants of the ad- 
vanced gravitation- wave detectors scheme which we op- 
timise, the model of technical noise and quantum noise 
which we use, and the optimisation procedure. In Sec. Ill, 
we discuss the results of optimisation. 



TABLE 1. The main parameters and their numerical values. 



xameter Value Description 
ujp 27rc/(1.064/im) Optical pump frequency 
M 40 kg Mirrors mass 

L 4 km Interferometer arms length 

/o < 500 W Input optical power 

Ic 200 kW or 840 kW Circulating power in the arms 
/abs < 1 W Power absorption in the ITMs 

and in the beamsplitter. 
?7 0.9 Photodetectors quantum 

efficiency 

/ 50 m Filter cavity length 

Af lOppm Filter cavity losses per bounce 



II. OPTIMISED CONFIGURATIONS AND THE 
OPTIMISATION PROCEDURE 

A. Interferometer configurations 

The following configurations are considered in this pa- 
per (see Fig. 1): 

• Plain: the ordinary signal- and power-recycled in- 
terferometer, similar to the Advanced LIGO, with 
vacuum input (no squeezing). 

• Squeezed: the same as the above, but with squeezed 
light injection into the dark port. 

• Pre-filtering: the same as previous, but with 
frequency-dependent squeezing angle implemented 
by means of the additional input filter cavity. 

• Post-filtering: the back action evading configura- 
tion with the additional output filter cavity and 
the squeezed light injection into the dark port. 

The last two configurations, in a more sophisticated 
two-cavities form, were first proposed in paper [25]. Here 
we consider more simple cases with only one relatively 
short filter cavity. Estimates (see, e.g^ [46]) show that 
in the presence of the technical noise, single filter cavity 
configurations provide almost the same sensitivity as the 
two-cavity ones, while being much less expensive in im- 
plementation. It was shown also in the paper [46] that the 
sensitivity of these phase filtering schemes is better than 
that of the amplitude filtering configuration [27, 31, 47]. 
Currently, it is these two schemes that are considered as 
the most probable variants for upgrading the Advanced 
LIGO. 

For the main parameters of the considered schemes, 
we assume the values close to the ones planned for the 
Advanced LIGO and its upgrades, see Table I. We sup- 
pose, however, that the input power can be higher, up 
to 500 W, and for the circulating power, along with the 
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FIG. 1. The schemes of the advanced gravitation- wave detector being considered in this paper. Top left: the signal- and 
power-recycled interferometer similar to the one planned for Advanced LIGO [35] (referred to as "plain"). Top right: the 
configuration with frequency-independent squeezed light injection into the dark port (referred to as "squeezed"). Bottom left: 
the configuration with injection of frequency-dependent squeezed light created by means of the input filter cavity (referred to as 
"pre- filtering" ) . Bottom right: the back action evading configuration with the additional output filter cavity and the squeezed 
light injection into the dark port (referred to as "post-filtering"). 
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"canonical" value of 840 kW, we consider the reduced one 
equal to 200 kW. 



B. Figures of merit 

The common figure of merit for the gravitational wave 
detectors sensitivity is the signal-to-noise ratio integral 
(SNR) for some standard gravitational waves source: 



/9^(p,7ViTM, A^ETm) 



^ 5^(1];7Vitm,A^etm;p) 2^ ' 

(2) 

where l^min and l^max are the minimal and the maximal 
frequencies of the GW detector sensitivity band, which 
we assume to be equal to 27r x 5Hz and 27r x 5 kHz, 
respectively; h{Q) is the gravitational wave strain signal 
spectrum; 

S^{Q] A/'iTM, ^ETM; P) = Equant (^5 ^ITM, ^ETM; P) 

+ St,,{n; iViTM, iVETM) + (3) 

is the spectral density of the sum noise normalized 
as the equivalent gravitation wave strain variation; 
'Equant 5 ^ITM, ^etm; p) is the quautum noise spectral 
density which depends on the numbers A^itm and A^etm 
of the coating layers doublets of the input and the end 
mirrors, respectively, as well as on the vector p of other 
optical parameters of the interferometer (we specify this 
vector below); 5'coat(^5 ^itm, ^etm) is the mirrors coat- 
ing thermal noise spectral density, which also depends 
on A/'iTM and A/etm; and 5't^ch(^) spectral 
density of the rest of technical noise sources. 

For our optimisation, we will use two standard types 
of the gravitational waves sources, (see, e.g, [48]). The 
first one is the GW signal from the inspiral stage of the 
binary neutron stars (BNS) collisions with 

I.. /om2 ^^<27rx 1.5kHz, 

10, Q > 271 X 1.5 kHz , 

which accounts mostly for the low-frequency noise. The 
second one is the gravitational wave "bursts" , with 



Burst 



K 



Burst 



(5) 



which require more broadband sensitivity from the detec- 
tor. Here i^BNS and i^Burst are some factors, depending 
on the astrophysical parameters of the signal source but 
not on the observation frequency Vt and detector optical 
parameters. 

Our treatment of the interferometer quantum noise is 
based on the works [25, 46, 49]. The corresponding ex- 
plicit equations for the quantum noise spectral densities 
'S'quant cumbcrsomc to be shown here, especially 

for the filter cavities based configurations, and can be 
found in Appendix A. 



In our calculation of the coating Brownian noise spec- 
tral density S'coat? follow the paper [17], which takes 
into account additional effects of light interference in the 
coating layers as well as the photoelastic effect. This 
method gives more precise correct estimate for the coat- 
ing Brownian noise spectral density, which is smaller than 
the previous estimates [12] by 3% - 13%, depending on 
coating layers number. 

For the value of spectral density S'^g^h other kinds 

of technical noise, we rely on the GWINC software tool 
[50]. 



C. Optimisation procedure 

For each of the four configurations, of the two figures 
of merit described above, and of the two values of circu- 
lating power, Ic = 840 kW and Ic = 200 kW (16 variants 
total), we found maximums of the SNR integrals (2) and 
(4) over the parameter vector p and in the numbers of 
coating layers doublets of ITM, A/'iTM, and ETM, A'etm- 

This optimisation was performed in two steps. First, 
we maximised the SNRs in p for each of A^itm and A^etm, 
using Nelder-Mead simplex method [51], and getting thus 
the semi-optimised values of SNR ^^(A/itm, A'etm) as a 
function of the numbers of coating layers doublets. At 
the second step, we found the optimal values of A^itm 
and A^ETM by a simple grid search. 

The following interferometer parameters were included 
into the parameter space of the vector p (depending on 
the configuration): 

Plain: p = {i?sRM, (/>sr, (/>lo}, (6a) 
Squeezed: p = {i?sRM, (/>SR, </>LO, e'', A}, (6b) 
Pre-/post-filtering: p = {i^sRM, ^SR, ^LO, e'', A, 7/, ^/}, 

(6c) 

where Rsrm is the signal recycling mirror power refiec- 
tivity, (j)sR is the signal recycling cavity detuning angle, 
^LO is the homodyne angle, is the squeezing factor 
and A is the squeezing angle, 7/ is the filter cavity half- 
bandwidth and Sf is the filter cavity detuning. We lim- 
ited the squeezing factor by 10 dB (e^ < vTO), according 
to the contemporary experimental achievements in the 
low- frequency squeezing [52, 53]. The input laser power 
was limited by 500 W. 

In Fig. 2, we draw a typical relative SNR gain factor 
as a function of ITM and ETM number of coating layers 
doublets defined as follows: 

r (AT AT \ - P^(^iTM, A'etm) 

^rel(AiTM, AetM j - T^nVdiT-^^TdiTT ' ^ ) 



where 



^^ITM — 



AT-def 
^^ETM 
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(8) 



are the default values for the number of coating layers 
doublets prescribed by GWINC for the Advanced LIGO. 
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FIG. 2. Relative gain in SNR as a function of the coating 
layers doublets numbers Nitm , A^etm • Cyan and magenta 
dots correspond to the default numbers A^itm? ^etm (8) given 
by GWINC. 



Cyan dots correspond to A^itm = ^vru^ 
genta ones — to A^etm = ^etm- ^^^y ^^^^ 
for each A^itm, the optimal value of A^etm exists, which 
corresponds to the balance of the quantum noise and 
the coating thermal noise. For smaller A/etm, the ETM 
transmittance is too big and the quantum noise domi- 
nates. Otherwise, the redundant number of the coating 
layers leads to a domination of the coating thermal noise. 

At the same time, the dependence of 
Grei(A/'iTM, ^etm) ou A/jtm displays no extremal 
behavior down to TVitm = 0. This means that our 
optimisation algorithm tries to remove the noisy input 
mirrors completely, switching to the GEO-600 pure 
Michelson configuration [6], and providing the necessary 
bandwidth of the interferometer by means of the signal 
recycling mirror. Unfortunately, in fact, this regime 
can not be considered as the optimal one due to the 
following reasons. First, in our consideration, we have 
not taken into account the noise introduced by the 
beamsplitter. In a standard Fabry-Perot /Michelson 
configuration with arm-cavity finesse T ^ 1, the 
beamsplitter-induced noise can be neglected, as each 
reflection of the pump laser light from the beamsplitter 
corresponds to T reflections inside the arm cavities, 
and thus it is suppressed by the factor J^~^. However, it 
becomes important if A^itm ^ and, therefore, ^ 1. 

In particular, in the Advanced LIGO, the motion of 
the beamsplitter relative to the ETMs will be monitored 
by means of an additional modulation sidebands of the 
carrier light resonant only in the power and signal re- 
cycling cavities, whereupon a sufficiently small value of 

will imply much more stringent requirements to this 
monitoring system [54]. 

Second, for a given value of the circulating optical 
power in the arm cavities /c, the power on the beamsplit- 
ter /bs is proportional to IcT~^, and very high power on 
the beamsplitter and ITMs can create undesirable effects 
like thermal lensing due to heating of the input mirrors 



and the beamsplitter by the absorbed optical power. 

Due to these reasons, we limited the minimal number 
of the coating layers doublets on ITM by A^itm = 4, 
which corresponds to ^ 20. This value is sufficiently 
high to neglect with good precision the noise of the beam- 
splitter and to provide a reasonably low optical power on 
the beamsplitter (in particular, in all our estimates, we 
limited the power absorbed in each of the input mirrors 
and in the beamsplitter by the value of IW). A bold 
red dot in Fig. 2 corresponds to the optimal configura- 
tion conditional on the above mentioned constraints. 

Another technical problem that may arise in the real 
interferometer is related to the pumping light leaving 
the arm cavities (that are usually impedance-matched for 
this light) through the ETMs. This light has a nonzero 
chance to return back into the cavities after scattering 
on the surrounding seismically not isolated objects and 
thus carrying random phase that may be a source of ad- 
ditional noise for the interferometer [55]. This scattering 
can be prevented by placing an absorbing plates with 
the power reflectivity i^abs behind the ETMs. Simple 
estimate shows that in the case of the specular (mirror- 
like) reflectivity of these plates, their positions have to 
be controlled with the precision about Tetm V Rahs ^ 
10~^-10~^ times relaxed compared with the signal dis- 
placement, where Tetm is the ETMs power transmissiv- 
ity (assuming i^abs ^ 0.01 and Tetm ~ lO^^-lO""^, see 
Tables IV, V). Note that in the double- mirror topology 
of [41], this factor is equal to just Tetm, that is, bigger 
by at least one order of magnitude. 



III. DISCUSSION 

The results of the optimisation are given in Tables II 
and III for the BNS and burst sources, respectively. In 
addition, the explicit values of the optimised vectors p 
are shown in Tables IV and V in the Appendix. 

It follows from these results, that reducing A^itm from 
A^i% = 8 to 4, and TVetm from N^^y^ = 19 to 14 for 
Ic = 840 kW and to 11-12 for 200 kW, it is possible to 
increase the signal to noise ratio for the BNS events 
by ~ 30%, and for the burst events — by 20%, see 
the corresponding columns labeled "Grei"- A few per 
cent surplus of Grei for low circulating power case {Ic = 
200 kW) over the high power one {Ic = 840 kW) is simply 
the result of a smaller number of coating layers required 
in the former case {N-^tm = 11-12 vs 14), i.e. the less 
circulating power is required for the same maximal input 
power, the more transparent can be the end mirrors. 

We are also able to compare the performance of the 
interferometers which employ different advanced tech- 
niques like squeezing, pre- and post-filtering with the 
baseline setup, characterized by 

/c = 840kW, r = 0, iVi%=8, Ni^^ = 19 . (9) 

For this purpose, we introduce the absolute gain defined 
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TABLE II. BNS optimisation. 



Configuration 


Ic 


A^ITM 


A^ETM 


/o [kW] 


/bs [kW] 






G'abs 


Plain 


200 kW 


4 


11 


0.5 


16 


- 


1.31 


1.24 


Squeezed 


200 kW 


4 


11 


0.5 


16 


1.5 


1.31 


1.25 


Pre-filtering 


200 kW 


4 


11 


0.5 


16 


10 


1.33 


2.13 


Post-filtering 


200 kW 


4 


12 


0.25 


16 


8.7 


1.33 


2.03 


Plain 


840 kW 


4 


14 


0.4 


68 






1.26 


Squeezed 


840 kW 


4 


14 


0.4 


68 


1.6 


1.26 


1.26 


Pre-filtering 


840 kW 


4 


14 


0.4 


68 


10 


1.30 


2.30 


Post-filtering 


840 kW 


4 


14 


0.4 


68 


10 


1.30 


2.25 



TABLE III. Bursts optimisation. 



Configuration 


Ic 


A^ITM 


A^ETM 


/o [kW] 


/bs [kW] 






G'abs 


Plain 


200 kW 


4 


11 


0.25 


16 




1.24 


0.59 


Squeezed 


200 kW 


4 


12 


0.25 


16 


10 


1.16 


1.22 


Pre-filtering 


200 kW 


4 


12 


0.25 


16 


10 


1.21 


1.41 


Post-filtering 


200 kW 


4 


12 


0.25 


16 


10 


1.21 


1.43 


Plain 


840 kW 


4 


14 


0.4 


68 






1.20 


Squeezed 


840 kW 


4 


14 


0.4 


68 


10 


1.17 


2.58 


Pre-filtering 


840 kW 


4 


14 


0.4 


68 


10 


1.21 


2.91 


Post-filtering 


840 kW 


4 


14 


0.4 


68 


10 


1.21 


2.92 



as: 

Gabs = , (10) 

^base 

with Pbase i"6lcited to the baseline interferometer. The 
corresponding values of Gabs are given in Tables II and 
III. 

These numbers tell us that for the GW bursts, the 
near-four-fold increase of the circulating power results 
in the two-fold relative rise of the SNR, yet almost the 
same result can be achieved by employing a 10 dB input 
squeezing. At the same time, the BNS signal-to noise 
ratio pIns i"emains almost unchanged with the power in- 
crease, and does not benefit noticeably from the squeez- 
ing (note that the optimisation algorithm virtually re- 
jects the squeezing, using only small value of r). 

The physics behind such a behavior is transparent. For 
the sensitivity of the detector to GW bursts depend on 
the high frequency shot-noise-dominated part of the noise 
budget (inversely proportional to the product of circulat- 
ing power Ic and squeezing factor e^^) to a greater ex- 
tent (/?'burst(^) 0^ Vt~^^'^) compared to the sensitivity to 
gravitational waves emitted by compact binary systems 
(^BNs(^) oc the former one displays the expected 

trend described above. At the same time, BNS sensitivity 
assigns more weight to the radiation-pressure-dominated 
low-frequency region and to medium frequencies equally 
sensitive to both parts of quantum noise, thus making it 
much harder to balance their influence, for the radiation 



pressure fluctuations have opposite dependence on power 
and squeezing factor as the shot noise. 

Another apparent conclusion yielding from our opti- 
misation is that even a relatively short (/ = 50 m) filter 
cavity is capable of quite significant sensitivity gain for 
both considered GW sources. Though there is no signifi- 
cant difference in performance for pre- and post-filtering 
schemes, given the assumed values of cavity length and 
losses per bounce, the pre-filtering scheme experimental 
implementation is easier and thus is more favorable. 

Plots of the optimal spectral densities for both quan- 
tum and technical noise are drawn in Figs. 3 and 4 for 
BNS and GW bursts, respectively. We elected for plot- 
ting only the most promising variants with fixed and 
frequency-dependent input squeezing. For comparison, 
the corresponding spectral densities for the default num- 
bers of the coating layers doublets (8), as well as for 
the "baseline" configuration (9) are shown there as well. 
These plots demonstrate that the optimisation procedure 
considered here, while suppressing noticeably the sum 
technical noise, introduces only minor changes into the 
quantum noise. The only exception is the BNS optimi- 
sation in the absence of the filter cavity, which reduces 
the quantum noise in the best sensitivity area around 
100 Hz, following the technical noise, at the price of nar- 
rower width of this region. 

We would like to emphasize that our solution to the 
coating thermal noise problem provides a 20-30% gain in 
the signal-to-noise ratio for the two most probable astro- 
physical sources of GWs which amounts in the 30-50% 
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FIG. 3. Optimal quantum and teclinical noise spectral densities for the BNS sources. Left column: schemes with frequency- 
independent squeezing (referred to as "squeezed"); right column: schemes with frequency-dependent input squeezing (referred 
to as "pre- filtering" ) . Top row: Ic = 200 kW, bottom row: Ic = 840 kW. Solid lines: optimised numbers of the coating 
layer doublets A^itm and A^etm; dashed lines: default numbers A^itm = 8 and A^etm = 19- Dash-dotted line: "baseline" 
configuration: Ic — 840 kW, A^itm = 8 and A^etm = 19, no squeezing; dotted line: SQL. 



increase of the event rate without any significant change 
to the optical layout of the GW interferometer. The price 
one has to pay, namely about two-fold increase of the in- 
put laser power, slightly higher absorbed power in the 
core optics (~ 1 W), and other issues discussed above, 
does not look too high in view of the potential benefits. 
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Appendix A: Quantum noise spectral densities 

1. Notations 

We base our consideration here on the Caves- 
Schumaker's two-photon formalism [56, 57]. Treatment 
of all of the considered schemes of the interferometers is 
done is accordance with [58] where all the detailed deriva- 
tions can be found. Two-photon quadrature vectors are 
denoted by boldface letters, and their components — by 
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FIG. 4. Optimal quantum and teclinical noise spectral densities 
independent squeezing (referred to as "squeezed"); right column: 
to as "pre- filtering" ) . Top row: Ic = 200 kW, bottom row: Ic 
layer doublets A^itm and A^etm; dashed lines: default numbers 
configuration, see Eq. (9); dotted line: SQL. 
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for the burst sources. Left column: schemes with frequency- 
schemes with frequency-dependent input squeezing (referred 
= 840 kW. Solid lines: optimised numbers of the coating 
NiTM — 8 and A^etm = 19- Dash-dotted line: "baseline" 



upright letters, e.g. 




2x2 matrices are denoted by "blackboard bold" letters, 
e.g. R. 

The quadrature amplitudes are normalized in such a 
way, that in the vacuum quantum state, their single-sided 
spectral densities are equal to one, e.g. 

5[gi] = 5[g2] = l. (A2) 

We use the following notation for the norms of the 
two-components quadrature vectors: 

VA=(^^^): ||A||2=aTa. (A3) 



2. Input-output relations 

a. Michelson/ Fabry- Perot interferometer 

Following the treatment of paper [49] (often referred to 
as "scaling law"), the Fourier- domain input/output rela- 
tions of a signal-recycled Michelson/Fabry-Perot inter- 
ferometer detector can be written in the equivalent form 
of the single effective Fabry-Perot cavity input /output 
relations: 

h{^) =R(l])a(l]) + T(l])g(l])+R(l])x(l^)^^. (A4) 

Here a, b, g are the quadrature vectors of, respectively, 
the input light, the output light and the vacuum noise 
which arises due to the optical losses in the interferome- 
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ter, 



i?ii R 



12 



^^"^ " v{vt) - js/n^ \^i?2i R22J ' 
= R22 = 271(7 - i^) - + js/n^ , 

R12 = -271 <5 , R21 = 271^ - 2J7i/f)2 , 
2\/7i72 / 7-if2 -<5 \ 



(A5a) 



(A5b) 



(A6) 



with i, 6, and q being the two-photon quadratures vec- 
tors of, respectively, the input, output and additional 
vacuum, induced by the optical losses in the filter cav- 
ity, quantum fields that fully describe all the filter cavity 
inputs and outputs, and 



Vf{n) ^Rf^, Rf^^J 

11 = Rf22 = 7/1 - 7/2 2injf2 

Rfl2 = —Rf2i = — 27/i(^/, 
2^7/17/2 f-ff-m -5f \ 



(A13a) 



(A13b) 



S = 



71 



71TM 



2 Vi^SRM sin(20SR)7lTM 
1 + 2 Vi^SRM COs(2(/)sr) + RsRM ' 

(1 — i^SRM) 7lTM 
1 + 2Vi^sRM COs(20sr) + i^SRM ' 

7 = 71 + 72 , 
c TiTM cTetm 



4L 



72 = 7ETM 



4L 



(A7a) 

(A7b) 
(A7c) 
(A7d) 



TiTM is the power transmissivity of the ITMs, Tetm is the 
ETMs transmissivity accumulating all the optical losses 
in the interferometer, i^sRM = 1 — ^srm is the power 
reflectivity of SRM, and (j)sR = cj^/sr/c is the single trip 
detuning phase of the carrier light in the SR-cavity with 
length IsR. 
Vector 



R(n) 



271 JM 1 / -S ' 



(A8) 



stands for the interferometer optical response function 
to the differential mechanical motion of the test masses 
(dARM mode) induced by the GW tidal force: 



G{n) = -MLn^h{n) 

acting on the ETMs, and 



(A9) 



(AlO) 



is the mechanical susceptibihty of the dARM mode mod- 
ified by the optical rigidity 



MJ5 



(All) 



b. Filter cavity 



Vf{n) = {jf-in)^^6}, (A14) 

^/ is the filter cavity detuning, 7j = 7^1 +7/2 is its 
half-bandwidth, 



f 



cAf 



(A15) 



Tf is the power transmissivity of the input / output mirror 
of filter cavity with length Lj, and ^/ is the coefficient 
of optical power losses per bounce. 



c. Homodyne detector 

We model the homodyne detector quantum efficiency 
7^ < 1 by an imaginary gray filter with a power transmis- 
sivity rj'. 



d' = V^d+ V'l-r/n, 



(A16) 



where d is the quadrature vector of the photodetector 
incident field, d' is the effective incoming field, and n 
is the additional vacuum noise associated with the the 
photodetector quantum inefficiency 1 — 77. 

The output signal of the homodyne detector (the pho- 
tocurrent) is proportional to 



^(^) oc H"^d'(n) oc H"^ [d(p) + edn(O)] , (A17) 



where 



and 



= a/ - - 1 



(A18) 



The input /output relations for the filter cavity have 
the following form: 



H 



cos 0LO 

sin 0LO 



(A19) 



6{Q) = M/(ri)i(ri) -h Tf{Q)q(Q), (A12) is the homodyne vector. 
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3. Configurations 



c. Pre-filtering 



In this case, 



a. "Plain^^ interferometer 



d = b, 



(A20) 



and the incoming field a is in the vacuum quantum state. 

Combination of Eqs. (A4, A20, A16) gives the following 
single-sided spectral density of the quantum noise in the 
units of gravitational- wave strain h: 



8 



'plain (^) M^L^n^ 

||Rt([])H| 



'^\\T\mr±4^ (A21) 



b. Squeezed input 

In this case, the incident field a is the result of squeez- 
ing of some vacuum filed z: 



with squeezing matrix defined as: 



(A22) 



g _ I cosh r + cos 2A sinh r sin 2A sinh r \ 
I sin 2A sinh r cosh r — cos 2A sinh r I 

(A23) 

Combining Eq. (A22) with Eqs. (A4, A20, A16), we ob- 
tain: 



lis: 



HII 



||Tt(f2)Hf + e2 



\H^R{n)xm\ 



. (A24) 



In the pre-filtering scheme, the incident field with 
frequency-dependent squeezing angle A(l^) is created by 
means of passing the squeezed light with frequency- 
independent squeezing (A22) through the filter cavity, 
see Eq. (A12): 

a{n) = Rf{n)Sz{n) + Tf{n)ci{n). (A25) 

Combining Eq. (A25) with (A4, A20, A16), we obtain: 



4( 



ISR^f (l^)Rt(f])Hf + ||Tj.(l^)Rt(n)H||2 



(A26) 



d. Post- filtering 



In the post-filtering scheme, the outgoing field of the 
interferometer passes though the filter cavity: 



(A27) 



The chain of Eqs. (A22, A4, A27, A16) gives the following 
spectral density: 



^post 



1 



\\T\{n)u\ 



2 , 2 



(A28) 
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TABLE IV. BNS optimisation. 



Configuration 


Ic 


TiTM 




i^SRM 


(/)SRC [rad] 


(/)LO [rad] 


A [rad] 




7/ [Hz] 


6f [Hz] 


Plain 


200kW 


0.15 


1 • 10"^ 


0.78 


0.30 


2.64 










Squeezed 
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0.15 


1 • 10~^ 


0.70 


0.29 


2.61 


-0.01 


1.5 


— 


— 


Pre-filtering 
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0.15 


1 • 10"^ 


0.07 


0.01 


1.29 


-0.33 


10 


158 


213 


Post-filtering 
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0.10 
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0.68 
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Plain 
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0.12 
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0.02 
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10 
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292 


Post-filtering 
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0.15 


1 • 10-^ 


0.06 
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1.77 


0.83 


10 


190 


264 
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0SRC [rad] 
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0.15 


1 • 10"^ 


0.67 


0.29 


2.06 










Squeezed 
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0.15 


6 • 10"^ 


0.01 


0.14 


1.60 


-0.01 


10 






Pre-filtering 
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6 • 10-^ 


0.04 


0.00 
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-0.06 


10 
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Post-filtering 
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0.15 


6 • 10-^ 


0.03 


0.00 
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0.04 


10 
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0.15 


1 • 10-^ 


0.38 


0.53 
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Squeezed 
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0.15 


1 • 10-^ 


0.05 
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-0.01 


10 






Pre-filtering 


840kW 
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10 


253 


254 
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0.03 


10 


260 
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